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primarily by low-energy spinwaves associated with cSDW order. We find that this 
inverse Fermi surface nesting mechanism is robust with respect to the addition of 
inter-orbital hopping. 
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1. Introduction 

The surprising discovery of high-temperature superconductivity in iron-pnictide 
compounds is perhaps the newest unsolved puzzle in condensed matter physics[T]. 
Calculations of the electronic band structure that include all five iron 3d orbitals, but 
that assume only weak inter-electron repulsion, predict two-dimensional (2D) nested 
Fermi surface pockets centered at zero momentum and at wavenumbers {TT/a)x and 
(7r/a)y[2]. This prediction is consistent with determinations of the electronic structure 
in iron-pnictide systems by angle-resolved photoemmission spectroscopy (ARPES) 
[3] |1] [5] [6] . Here a denotes the lattice constant of the square lattice of iron atoms that 
stacks up to form iron-pnictide systems. Band structure calculations also predict an 
ordered magnetic moment of 2 Bohr magnetons (hb) in the commensurate spin-density 
wave (cSDW) state that exists in parent compounds at low temperature. The ordered 
cSDW moment that is measured in iron-pnictide systems by elastic neutron diffraction 
can be much smaller [7]. It reaches values as low as 0.3 /is . Frustrated Heisenberg 
models that assume local magnetic moments at each iron atom can successfully 
account for the weak cSDW that exists in parent compounds [8] [9] [TO] |1 1 j . on the other 
hand. They can also can give a good account of the low-energy spin excitations 
near cSDW momenta that have been uncovered in iron-pnictide systems by inelastic 
neutron scattering [12] fi3\ fi5\ flQ\ [T7] . Such Heisenberg models have an insulating 
groundstate, however, that is a result of strong inter-electron repulsion. This fact 
conflicts with the metallic nature of iron-pnictide superconductors and their parent 
compounds. 

Iron-pnictide high-temperature superconductors then appear to be simultaneously 
itinerant-electron metals and local-moment magnets. We identify a way to resolve 
this dilemma from the limit of strong inter-electron repulsion by injecting a low 
concentration of mobile charge carriers into a local-moment cSDW that experiences 
magnetic frustration that is off-diagonal in the orbital indices|9]. Both a Schwinger- 
boson-slave-fermion mean-field theory and exact computer calculations on a finite-size 
system find a quantum phase transition into a hidden half metal state as Hund's 
Rule coupling weakens [18]. The new groundstate shows long range antiferromagnetic 
order across the d+ = 3d(^x+iy)z and d— = 3d(x-iy)z orbitals. We use the d± orbitals 
because they are the least localized ones, and because they hence maximize Hund's Rule 
coupling. This orbital basis therefore maximizes the tendency to form local magnetic 
moments on each iron atom. The new hidden half metal state exhibits coherent intra- 
orbital hole motion, which yields two Fermi surface hole pockets centered at zero 2D 
momentum. Inter-orbital hopping fails to destroy the state. Instead, it can lead to 
an orbitally sensitive reverse Fermi-surface nesting mechanism near the quantum phase 
transition that is mediated by zero-energy cSDW spin-wave excitations. (See Figs. [9] 
and [121) The latter, in particular, can result in a nested "electron" Fermi surface pocket 
at wavenumber {tt /a)x{y) with predominantly 3dyz{3dxz) orbital character. Mean-field 
theory predicts that the single-particle excitations here show a mixture of electron and 
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hole dispersion |1 8]. We now pass to the analysis of the theoretical model. 
2. Two-Orbital t-J Model 

We shall now introduce a two-orbital t-J model that describes the low-energy electronic 
excitations in iron-pnictide high-temperature superconductors [18]. Spin-1/2 moments 
exist over the 3(i(j.+jj^)^ and 3d(^x-iy)z orbitals at each iron atom. These are the least 
localized orbitals within the 2D subspace spanned by the degenerate Sd^z and Sdy^ 
orbitals, and they therefore result in the largest possible Hund's Rule coupling constant, 
—Jo- This fact is demonstrated explicitly in the Appendix for the case of hydrogenic 
3d orbitals with the bare Coulomb interaction. In contrast, the commonly used 3dxz 
and 3dyz orbitals from Chemistry are the most localized ones within the same 2D 
subspace [19] |20] [21], and they result in the smallest possible Hund's Rule coupling. 
The isotropic nature of the 3d(x±iy)z orbitals also implies isotropic Heisenberg exchange 
coupling constants Ji and J2 across nearest-neighbor and next-nearest-neighbor links 
of the square lattice of iron atoms within each iron-pnictide layer. The two-orbital t-J 
model over a square-lattice of iron atoms in an isolated layer then reads 

^ = - E E T.('T'ci,sc,,, + h.c.) J2it-2''i./c,,, + h.c.) 

+ 2"^o E [E "^^'^ + E E -^I'^^i'^ ' ^i^P + E E -^i'^^ha ■ 

i a a,0 a,l3 

(1) 

Above, Si^a is the spin operator that acts on the spin sq = 1/2 state localized on orbital 
a at site i, while 4,a,s is the corresponding electron creation operator that takes into 
account the constraint against double-occupancy at a site-orbital. The orbitals a are 
either d+ = 3d(^x+iy)z or d— = 3d(^x_iy)z , while the sites i run over the square lattice 
of iron atoms. The isotropy and the degeneracy of this pair of orbital states yields two 
independent and isotropic Heisenberg exchange coupling constants for nearest-neighbor 
and for next-nearest-neighbor links, {i,j) and {{i,j)), respectively, n = 1 and 2: 

= j\\ = and J:^+'^- = J^ = J:^-'^\ (2) 

Correlated hopping of an electron in orbital /3 to a neighboring unoccupied orbital a 
is controlled by the hopping matrix elements t^'^ and ^2'^^. Let us take real hopping 
matrix elements in the {3drcz,3dyz) orbital basis |22]. The hopping matrix elements in 
the present d± basis are then given by 

C''^ = I [tr + 1^] ± I K'' - tr] , (3 a) 

The symmetry relation t^'^^-^ = t^g) in the {3dxz, 'idy^) orbital basis then yields diagonal 
hopping matrix elements that are real and isotropic in the {d+, d—) basis: 

tf = tf and lmtf'''± = 0, (4a) 
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tf^'^^ix + y) =t'!t' {y-x) and Imtf '''^ = 0. (46) 

The off-diagonal hopping matrix elements tf^''^^ have d-wave symmetry, on the other 
hand. Also, the identities ti'^ = = t\'^ yields real off-diagonal hopping matrix elements 
across nearest neighbors, while the symmetry relation ^2'^ = ^2'^ yields pure-imaginary 
off-diagonal hopping matrix elements across next-nearest neighbors: 

tf ''^^(a;) = -tf '^^(^/) and Imtf ''^^ = 0, (5a) 

tf '''^(a; + J/) = -tf '^^(jZ-i) and Retf '^^ = 0. {5b) 

Henceforth, we shall change the notation for the hopping matrix elements to t^^^ = 

.d±,d± . ,_L _ ,d±,dT J f _,d-,d+ _ j.± _ ^d+4- 
1(2) ' 1 — 1 ' aiiu. 10 — 62 — ^2 — ^2 

Last, it is useful to point out that because t^ is real, a global swap of the orbitals, 
(i± — i- (i=i=, is an exact symmetry of the two-orbital t-J model ([T]) in the absence of 
next-nearest neighbor inter-orbital hopping: = 0. Let P^^j denote the global swap 
operation of the d+ and d— orbitals. Eigenstates of the t-J model Hamiltonian ([1]) are 
then even (+) or odd (— ) under it, with respective forms ± P^j|\l/). Further, in the 
limit Jo — i- —00 where Hund's Rule is obeyed, a spin triplet state exists at iron atoms 
that do not have any holes. Such spin-1 local moments are clearly even under P^^j. This 
implies that a one-hole state that is even under P^g has a local Sd^^ orbital when Hund's 
Rule is obeyed, and that a one-hole state that is odd under P^j has a local ?idyz orbital! 
At the opposite extreme Jq — )■ +00 where Hund's Rule is maximally violated, a spin 
singlet state exists at iron atoms without any holes. Such spin-0 iron sites are clearly 
odd under P^j. In the case of an even number of iron sites, this implies that a one-hole 
state that is even under P^j instead has a local "idy^ orbital, and that a one-hole state 
that is odd under P^j instead has a local Sc/^.^ orbital! One-hole states that are even 
or odd under P^^j must then have a mixture of ?)dxz and ?)dyz orbital character at finite 
Hund's Rule coupling by continuity. 

And because is pure imaginary, the global orbital swap operation is an exact 
symmetry in the absence of nearest neighbor inter-orbital hopping, t^ = 0, after making 
the global gauge transformation of the orbitals \d±) i^^^'^\d±). Eigenstates of the t-J 
model Hamiltonian ([1]) are then even (+) or odd (— ) under the combined operation P^j, 
where they have the form ± P^j|\l/). Notice that the global gauge transformation 
rotates the orbital axes by 45 degrees. Repeating the previous arguments yields that 
a one-hole state that is even under P^j has a local 3dx'z orbital when Hund's Rule is 
obeyed, and that a one-hole state that is odd under P^j has a local 3dy'z orbital in such 
case. Here {x', y') denote the planar orbital coordinates along the next-nearest neighbor 
links. 

3. Hidden Half Metal State 

In general, the Heisenberg exchange constants have a direct ferromagnetic contribution 
from the exchange Coulomb integral (see Appendix) and an indirect antiferromagnetic 
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contribution from the super-exchange mechanism through the pnictide atom [23]: 

■Ji{2) = -Jm (drct) + Jf(f) (sprx) . (6) 
We shall assume that the super-exchange contribution is independent of the orbital 
indices: J^^^(sprx) = J^f^^^'' ■ This is the case if the pnictide orbital in question is the 
Px, Py, or pz one, for example. Recent density-functional calculations for the electronic 
structure of iron-pnictide materials find that direct exchange and super exchange can 
cancel out across nearest neighbors |24j |25j . We shall assume that this is the case for 
nearest-neighbor Heisenberg exchange across the same orbital: j| = 0. It is further 
reasonable to assume that J^2)'^^(drct) is smaller in magnitude than J^2j'^^(drct). This 
assumption is borne out in the limit of large overlap between neighboring iron 3d orbitals 
(see Appendix), where —jf^^'^^ (diet) is 10 times smaller than —jf^^^'^^dicty. We are 
then left with a net antiferromagnetic Heisenberg exchange across nearest neighbors 
between different orbitals due to super-exchange: > 0. Last, the previous DFT 
calculations find that J^''^(drct) is negligible. This leaves antiferromagnetic Heisenberg 
exchange across next-nearest neighbors due to super-exchange, which we assume is 
invariant with respect to a rotation of the orbital basis: j| = > 0. 

In the absence of mobile holes, the off-diagonal magnetic frustration described by 
the above list of Heisenberg exchange constants results in a commensurate spin-density 
wave (cSDW) along one of the principal axes of the square lattice of iron atoms at 
strong Hund's Rule coupling, —Jo, and for J2 > Ji/2. Here Ji(2) = (jf(2) + '^i(2))/2- 
An antiferromagnetic state with opposing ferromagnetic order over the d-\- and the d— 
orbitals appears at Hund's Rule coupling below a critical value of —Jqc = 2( Jf*" — 4j| 
in the large-So limit [11]. The latter marks a quantum-critical phase transition between 
the cSDW at strong Hund's Rule coupling and the hidden ferromagnet at weak Hund's 
Rule coupling. (Cf. Fig. [H) Observe now that injecting mobile holes while inter-orbital 
hopping is suppressed preserves the \d+\d- spin order of the hidden ferromagnetic 
state in the classical limit [TT]. This is a semi-classical picture of a hidden half metal 
state with coherent intra-orbital propagation of holes that follows the dispersion relation 

Eeik) = — 2t!|(cos kj-a + cos kya) — 2t|(cos A;+a + cos k_a), (7) 

where k± = kx ± ky. It implies two degenerate hole Fermi surface pockets centered at 
zero momentum when t\ < and tf + 2t| < 0, in qualitative agreement with ARPES on 
iron-pnictide superconductors. Below, we shall demonstrate that the hidden half-metal 
state is robust with respect to the addition of inter-orbital hopping. We will also show 
that a reverse SDW nesting effect exists at the quantum-critical point that separates 
the hidden half-metal state from the cSDW state. It results in copies of the hole Fermi- 
surface pockets centered at cSDW wavenumbers {n/a)x and {7i/a)y. We shall see that 
the latter have mixed electron and hole character according to meanfield theory. 

3.1. Schwinger-Boson-Slave-Fermion Mean- field Theory 

Unlike the doped Neel state in the one-orbital nearest-neighbor t-J model for copper- 
oxygen planes in high-Tc superconductors [26] [27] . hole propagation in the hidden half 
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Figure 1. Shown is tlie phase diagram obtained from the Schwinger-boson-slave- 
fermion mean- field theory of the two-orbital t-J model, Eq. ([T]). The red line marks 
a quantum-critical point that separates a charge-carrier poor cSDW metal from a 
charge-carrier rich hidden half metal at fixed Hund's Rule coupling. 



metal does not produce strings of overturned spins that disrupt the antiferromagnetic 
order when inter-orbital hopping is absent. This suggests that the mean- field 
approximation for the Schwinger-boson-slave-fermion formulation of the two-orbital t- 
J model (II]) is applicable for this state at weak inter-orbital hopping. The constraint 
against double occupancy is enforced by expressing the correlated electron annihilation 
operator as Ci^a,s = bi^a,sfia along with the new constraint 

per site i, per orbital a, where sq = 1/2 is the electron spin |28][29]. Here, bi^a,t 
and bi^a,i are the annihilation operators for the pair of Schwinger bosons, while fi^a is 
the annihilation operator for the spinless slave fermion. The spin operator is then 
expressed in the usual way: Si^^ = l^Y.s,sf fi,o'bla,s^s,s'bi,a,s'fla- Henceforth, we 
shall average the dynamics of the slave fermions over the bulk by writing Si^a — 
(1 — x)^h^^ ^^bj ^ g<Ts^s'bi,a,s' instead. Here, x denotes the concentration of holes 
per orbital. Following Arovas and Auerbach|30] |31j . we introduce symmetric versus 
antisymmetric meanfields with respect to spin fiip for the ferromagnetic versus the 
antiferromagnetic links of the hidden ferromagnetic state (Fig. [1]): 

{Kd±A'i±^i) = Qi = (^U±,t^^-.'^±.t) in = 1, 2), (9) 
Mbi4±,ibj,d^,t) = Qn = T{bi4±,tbj^dT,i) = 0, 1, 2), (10) 

where j = i for on-site links n = 0, where j = i + x{y) for nearest-neighbor links n = 1, 
and where j = i ± x -\- y for next-nearest-neighbor links n = 2. Next, we introduce 
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l{fUfj4-r = Pt = lifU-hd^) (n = 1,2). 



(11) 

i,d+J3,a-/ n\Ji.d-JJ,a+/ V" (12) 

Last, we introduce meanfields associated with hopping of the Schwinger bosons across 
different orbitals: 

{bid+,sb,4-,sr = K = {bld-,sh^+,s) {n =1,2), (13) 
where s =t, i- Now assume that all of the mean fields are homogeneous, and assume that 
the mean fields -P]'(2)5 Qi(2) Qi{2) isotropic. The mean-field approximation for the 
t-J model Hamiltonian (P) then has the form Hmf = Ho[P, Q, R] + hI'^^ + ijf"^ + Hf, 
where 

-f^o = 2 ^ JqQqQq* + 

i 

2 J2 5Z[-(24q!^J' + c.c.) - {2tiRiPi + c.c.) - j'}Q\Qf + J^Q^Q^] + 

2 H E[-(2tIg|p| + c.c.) - {2t^R^P^ + c.c.) - 4q\q\* + J'^^Q^Qi*] 

{{hi)) " 

consolidates the bilinear terms among the mean fields, where 



(+) 



and 



(-) 





t - 






bj,+,i 
. bd-,i _ 




bd+,i 


t - 


bd-,i 
P 

"d+,t 





fill 



+fij 



fiii 





+fix 
-fi, 



fill 
fi'r 



fi'^ 

fill 



fill 







-fix 


fi'^ 


fill 


+fix 








+fi± 


fill 


fiY 


-fi. 





fi'r 


fill 

















. "d-,i. _ 




bd+,i 




bd~,i 




V 







are the pair of Hamiltonians for free Schwinger bosons, with matrix elements 
fi||(fe)= A + 4^(jM + 2tiP^)7n(fc), 

n=l,2 

fi^(fc) = ZnJ'n^Qnlnik), 
n=0,l,2 

^'±if^) =^Y1 ^^(^)^i^(^) cos(fc„a) 



(14) 
(15) 



n=x,y 



+ 4j2ti{y ± a;)P2^*(i' ± i) cos(fc±a). 



(16) 
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is the Hamiltonian for free slave fermions, with matrix elements 



(17) 



n=l,2 



e±{k) = 4 ^ t^{n)R^{n)cos{kna) 

n=x,y 

+ 4 ^4(j/ d= x)R^{y ± x) cos(fc±a). (18) 
± 

Above, the destruction operators for Schwinger bosons of momentum ±fe are defined by 
ba,s{k) = Npe~^^'^ ^■e'^''''^'bi^a,s and ba,s{k) = fea,s(— fc), while the destruction operator 



for a slave fermion of momentum k is fa{k) = A^Pe" Yli^^ ^' fi,a- Also above, 
zq = 1 and Zi(2) = 4 give the coordination number, and we define 7o(fe) = 1, 
7i(fc) = ^{cos k^a + cos kya), and 72(fc) = |(cosA;+a + cos/c_a), where A;± = fc^; ± /c^,. 
Last, A is the boson chemical potential that enforces the constraint against double 
occupancy (|8]) on average over the bulk of the system, and the effect of mobile holes 
on the Heisenberg spin-exchange is accounted for by the effective exchange coupling 
constants [29] J' = (1 — x)^ J. 

The spin-excitation spectrum of the hidden half metal state is obtained from the 
sum of the Schwinger bosons Hamiltonians, Hj^'^^ -\- Hj^'^ , in two steps. We first make a 
two-orbital Bogoliubov transformation of the boson field: 

ba,s = (cosh 9)(3a,s — (sgu a) (sgu s) (siuh 9)f3l,^s, (19a) 
%,s = (cosh 9)f3l^g - (sgna)(sgns)(sinh 9)f3a,s, (196) 

with cosh 29 = and sinh 29 = n^/ui°\ where u^^^ = (1]J - ^1^/^. Here, we 

use the notation d± = ci=F in the orbital index and s = — s in the spin index. The sum 
nj^'^^ -\- hI~^ then transforms to 



k s=t,4. 
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Let e^*'^*'-'^^ denote the phase factor of the inter-orbital matrix element n'j_(k), and let 
ko = 0, TT denote bonding (-I-) and anti-bonding (— ) superpositions among the d± 
orbitals after making the gauge transformation e^^^''^'^\ The similarity transform 

Mk) = e^^^^(^) [2-1/2/3,(0, k) T T^'^Ut^. k)\ 
then reduces the Schwinger boson Hamiltonian to 

= ^ E E E ^Amimsi}^) + fis{-km 



(20) 
(21) 



with the spectrum 

cU{,(A;o, k) 



(22) 
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The charge-excitation spectrum due to the slave fermions, on the other hand, is obtained 
directly by the similarity transform 

= e±*^/(^)[2-V2/(0,fc) T2-V2/(7r,fc)], (23) 

where e'^^^f^''^ denotes the phase factor of the inter-orbital matrix element e±{k). It 
yields the diagonal form 

E J2'fik)fHk)f{k) (24) 

fco=0,7r fc 

for the slave-fermion Hamiltonian, with spectrum 

efiko,k)=enik) + e'^^'\e^ik)\. (25) 

The low-energy spin and charge excitations that result from Ub and Ef will be discussed 
in detail below. 

We must first obtain the mean fields P, Q and R, however. Taking the quantum- 
thermal average of the constraint against double occupancy (IHl), and averaging it over 
the bulk, we get the principal mean-field equation 

So + ^{l-x)= Af-' 5^(cosh 26) (1 + nB[cob{k)]), (26) 

k 

where = 2A''Fe, and where ub denotes the Bose-Einstein distribution: UBiu)) = 
[exp{u/kBT) — 1]^^. Next, averaging the definitions of the mean fields (!9|)- (fT3|) over the 
bulk yields self-consistent equations 

Qi = A/""' $^7n(fc)[cosh 2e{k)] + UBluJbik)]) (n = 1, 2) (27a) 

k 

Qn = A/""' Yl ^"(^) [^^^^ + M^bik)]) (n = 0, 1, 2) (276) 



and 



Ri{h) = Af-^ J](cos kna)e''"'e^'^''^''^nB[ub{k)] (n = x,y) 

k 

R^{y ±x)= AT-i ^{cos k±a)e'^'e^'^''^^^nB[uJi,{k)] (28) 



k 



associated with the Schwinger bosons, and it yields self-consistent equations 

Pi = l^f-'J2^n{kW[ef{k)] (n= 1,2) (29) 

k 

and 

P^^(n) = -A/"-! J2^cos kna)e"''>e^'^f^'''^nF[€f{k)] {n = x, y) 

k 

P^{y ±x) = -U-^ ^(cos k±a)e'^°e^'^f^^^nF[e f{k)] (30) 

k 

associated with the slave fermions. Above, up denotes the Fermi-Dirac distribution: 
np^e) = (exp[(£ — ^)/kBT] -\- 1)~^, with a chemical potential fi. Let us now assume that 
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the mean fields Rf- and i?^ are both isotropic, which will be shown to be self consistent. 
This yields the form 

£±{k) = At^{x)R^[cos{kxa) — cos(fcya)] 

+ 4t^(i + ^)i?^[cos(A;+a) - cos(A;_a)], (31) 

for the inter-orbital matrix element experienced by slave fermions, where tj^ is real and 
where is pure imaginary. It alternates in sign when k is rotated by 7r/2. This d-wave 
symmetry implies that both P^{x) -\- P-^{y) = J^^^ ^j^e^''°e'^^^f^'^'''yi{k)nF[e f{k)] and 
P^{x + y) + P^{y -x) = J\f-'^ Y.k e'''"e2%(*')72(fc)ni7[£:/(A;)] vanish. Their mean-field 
equations then reduce to 

P^(x) = iAr-i^e''=«e2^^^('=)[(cos Ka) - (cos kya)]nF[ef{k)], 

k 

P^{x + y) = ^A/"-^ J]e*'=«e2*'^^^'^)[(cos k+a) - (cos k_a)]nF[ef{k)], (32) 

k 

with opposite signs for the orthogonal links. The latter d-wave symmetry then yields 
the form 

fl'^ik) = 8ti{x)P,^*{x)-f,{k) + Stiix + y)P^*{x + t/)72(fc) (33) 

for the inter-orbital matrix element experienced by Schwinger bosons. It has s-wave 
symmetry, and it is thus consistent with the previous assumptions. 

To proceed further, we shall assume that both inter-orbital mean fields for 
Schwinger bosons R^ and i?^ are real and positive. This will also be confirmed self- 
consistently. The slave-fermion matrix elements are then approximated by 

e\\{k) = e||(0) -t|||fc|V (34a) 
e±ik) = - [txcos(20) +it'^sin(20)]|fc|V (346) 

for k near zero, where = 2t|Ql + 4t|Q|, where t± = 2ti{x)R^, and where 
t'j_ = 4t^(A + y)R2 /i are real hopping amplitudes. Above, (p denotes the angle that 
k makes with x axis. Henceforth, we shall assume hole bands: ty < 0. The constant 
energy contours ep = £\\{k) ± |e_L(fc)| then yield Fermi surfaces 

4±(0)a2 = ep/{-t\\ ± [tl cos'(20) + t'l ?,in^{2(P)Y'^) (35) 

at low doping a; ^ 1, where e'p = ep — £||(0). These are shown in Fig. O The 
average of the volume of phase space contained by the inner (-I-) Fermi surface 
and the outer (— ) Fermi surface gives the hole concentration per orbital: x = 
(27r)^^ Jq"^ d(j) ^_|_ Jq'"'''^'^^ dk k a^/2. The resulting angle integral is obtained by applying 
the residue theorem after making the change of variable z = e^*"'^, which yields the 
relationship 

Attx = eU-t||)/(t| - tlY^'itl - t'ly/^ (36) 

between the Fermi energy and the concentration of mobile holes. Next, the intra-layer 
hopping fields for slave fermions are given by Pj = (27r)^^ f^"^ dcj) J2± /o^^'''^'^^ dkk{l — 
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Figure 2. Shown are the hole Fermi surfaces centered at zero 2D momentum. The 
concentration of holes is universally set to a; = 0.01, while the hopping matrix elements 
from left to right are (a) ti_ = 0.40(-t||) and t'j^ = 0.23(-t||) , (b) ti_ = 0.28(-t||) and 
4 = 0.28(-t||) , and (c) t± = 0.23(-t||) and t'j_ = 0.40(-t||) . Above, orbital labels 
are only approximate near the points on the Fermi surfaces at which they change. 



^/c^a^)a^/4 at X ^ 1. Performing the integral after the previous change of variable 
yields the result 

pii = 1 _ 2 ^ ~ ^i^±A5 ,o7N 

for these amplitudes. Last, the mean fields for inter-orbital hopping by slave fermions are 
given by P^{x) = (27r)"2 J^^ d(j) J^^^^'^^ dk ^ {t cos 20)e2%('^) aVl6 and by P^{x+ 
y) = {27t)-^ J^^ d(f)J2± lo'^^^'^^ dk k^Tsin 2(j))e^'^f^'''^ a^/8 at low hole concentrations, 
where 

„2%(fc) ^ t± cos(20) + it'^ sin(20) 

[ticos2(20) + t'_2sin2(20)]i/2 ^ ' 

is the phase shift at k near zero. Employing again the change of variable z = e^^"^ in the 
resulting angle integrals then yields the following expressions for these amplitudes: 



TT 



t\ ni-tw/' 
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Substitution into the expression for the inter-orbital matrix element experienced by 
Schwinger bosons fl33l) reduces it to 

|9.±|2 /t2 _ 1/2 I , _L|2 - ^2 N 1/2 

r]l(fc)=vra:X^(^) 7i(fc) + vrx^i?2^^ (/^^j l2{k). 

(41) 

Observe that the phase factor associated with Schwinger bosons is then e'^^^''^'^^ = 
sgnt|| = — 1 at near 0. The principal mean-field equation (!26|) implies ideal Bose- 
Einstein condensation (BEG) of the Schwinger bosons into the bottom of the spectrum 
f l22p at 3 momentum (/cq, k) = {it, 0) in the zero-temperature limit. Comparison of (12^ 
with the mean-field equations fl27fflP and (\27b\i yields the mean-field values Qn = Sq 
and Q;^ = Si, in the large-So limit, where Si = Sotanh 26^(0) = [fi_|_(0)/n||(0)]so. Ideal 
BEC necessarily requires Ub{7i,0) = 0. This yields ultimately that tanh 26* (0) = (1 + 
[fi'j^(0)/f2x(0)]^)~^/^. Inter-orbital hopping therefore diminishes the antiferromagnetic 
order in the hidden half metal state: Si < Sq. Notice, however, that the correction 
is small and of order x^! Last, ideal BEC implies a unique mean-field amplitude 
= S2 by meanfield equations where ss = so/cosh 2^(0) = [\Q'^{0)\/Q±{0)]/ {1 + 
[n'_j_(0)/f2_L(0)]^)^/^. Observe now that the latter ratio of frequencies takes the form 
|ri'^(0)|/fi_L(0) = Ao{s2/soSi) as X —J- 0, in which case si — )■ sq. Here, Aq is a positive 
constant that is small compared to unity at low hole concentration. Substituting this 
form into the previous definition of S2 yields only the trivial solution S2 = 0. We 
therefore conclude that inter-orbital hopping is dynamically suppressed in the present 
Schwinger-boson-slave-fermion mean-field theory for the hidden half-metal state at large 
electron spin sq and low hole concentration. In particular, Eqs. ( I39l) and ( l40ll imply 
that both the inter-orbital fields and i?^ are null. 

Dynamical suppression of inter-orbital hopping at large-so implies a null inter- 
orbital hopping matrix element for Schwinger bosons by Eq. ( HT]) : Q'j_{k) = 0. 
Inspection of the spectrum for Schwinger bosons (122|) then yields that it disperses as 

uJb{ko, k) = vo\k\ near k = 0, (42) 

where 

vo = 2soa(l - xf{[J^ - J| (x) + 2J^ - 2 jj (x)] ■ [^^o + + 2J^]y/\43) 

is the longitudinal spin-wave velocity, with Jn{x) = jj| + (1 — x)~'^SQ^tnX for n = 1,2. 
At k near cSDW wavenumbers (vr/a, 0) and (0, vr/a), the spectrum for Schwinger bosons 
(1221) disperses anisotropically as 

Ub{ko, k) = [AlsDw + - n/af + v^{kt/-fcSDwfY^\ (44) 
with a spin gap 

AcSDW = (1 - x)\2so)[{U^ - Joc)(Jo - Joc)f\ (45) 
and with an anisotropy parameter 

IcSDW = ([2j|(x) + 24 + jf(x) + Ji^]/[2j|(x) + 24-jjl(x)-Ji^])i/2,(46) 
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which is greater than unity. Here, the critical Hund's Rule coupling at which A^sdw 
is 

- Joc = 2{J^ -4)- 44' - (1 - x)-^SQ^2t\\x. (47) 

The hidden half-metal is stable at weak to moderate Hund's Rule coupling — Jq < ~Joc- 
Equation ( 147|) therefore implies that intra-orbital hole hopping stabilizes the hidden 
half- metal state for hole bands, < 0. (See Fig. [H) We conclude that the hidden 
half-metal is robust with respect to the presence of inter-orbital hopping within the 
present Schwinger-boson-slave-fermion mean-field theory. 

Finally, returning to the generic unoptimized mean-field theory, the inversion of the 
similarity transform ( 1231) yields slave-fermion states that are annihilated by 

+i5f{k) -i5f{k) 

f{ko, k) = -^^f,_(k) + e^^-^^fa^. (48) 

Their corresponding Wannier wave functions then depend on the azimuthal angle (j)' 
about each iron atom as cos[0' — 5/ (A;)] in the even channel and as sin[(/)' — 5/(fc)] in the 
odd channel. Figure |2] depicts slave-fermion Fermi surfaces fl35|) at low doping x = 0.01. 
Dynamical suppression of inter-orbital hopping implies that the matrix element for 
inter-orbital hopping of slave fermions ( !T8|) is null, however. The inner and outer Fermi 
surface hole pockets depicted by Fig. [2] must therefore collapse into two degenerate 
circular Fermi surface hole pockets at low doping, a; ^ 1, at large spin sq. We believe 
that this effect is due to the inability of the two-orbital t-J model ([1]) to "erase" strings 
of overturned spins that result from inter-orbital hole hopping in the classical limit, at 
large spin so[26][27j. As we will see in the next section, this limiting result persists for 
true electron spin sq = 1/2 in certain cases. 

3.2. Spinwaves 

We shall now determine the spin excitations of the hidden half metal by directly 
computing the dynamical spin-spin correlation function from the Schwinger-boson- 
slave-fermion mean-field theory. First, we identify the true (/cq = 0) versus the 
hidden (fco = tt) spin at each iron atom i: Si{ko) = Si^d- + G^'"^Si^d+- Consider 
then the transverse dynamical spin-spin correlation function (^(fco) ■ '^"(^o))!"'" = 
|(S'"'"(A;o)5"~(A;o) + S~{ko)S'~^{kQ)) at space-time points {r,t) and {r',t'). Averaging 
the slave-fermion dynamics over the bulk, we obtain the form 

(^(^o) • s'{ko))\t = (1 - E E ^^'"^""'^ ■ 

a=0,l ,9=0,1 

for its Fourier transform, where iG^^^^{ri,t;rj,t') = {hi,a,s{t)b] p s{t')) and 
iF^'j^{ri,t;rj,t') = (&j,a,s(^)&j,,g,s(^')) are the regular and the anomalous Greens func- 
tions for the Schwinger bosons, and where the notation f * g denotes a convolution in 
frequency and momentum. Above, we use the index and 1 for the orbitals d— and d-\-. 



Fermi- Surface Nesting in Reverse from Local-Moment Description 



14 




Figure 3. The mean-field result for the imaginary part of the dynamical spin response 
function in the zero-temperature-large-sg limit, Eqs. (|5ip and (|52p . is evaluated with 
the following set of parameters: j| = 0, Ji > 0, = 0.3 Jf^ — J^^, if — — 5Jf^, 
t\ = 0, X = 0.01, and sq = 1/2. Dynamical suppression of inter-orbital hopping is 
presumed. The Hund's Rule coupling is set to — Jq ~ —Joc — O.lJf^ just off the QCP 
in panel (A), while it is set to the critical value — Jqc in panel (B). Low-energy contours 
are displayed in the inset. 



Also, henceforth in this subsection, we shall identify ko = and vr with the respective 
labels -|- and — . After substitution of the Bogoliubov transformation fll9ffl[ &) and of 
the similarity transform fl20|) . a standard summation of the Matsubara frequencies in 
the convolution yields the following Auerbach-Arovas expression for the dynamical spin 
correlator at T > 



90=+,- q 



^{S{±)■S'{±))\i^ = ^^-^-Jf^ ^ X^|(l + cosh[2%)T2e(q-fc)]) 

B[uJb{qo, q)] - nB[ujb{qo±, q - 

V - Ubiqo, q) + UJb{qQ±, q-k 



E nsluJbiqo, q)] - nB[uJb{qo±, q - k)] ^ 
;y — Ml ( -i- (,), ( nr.-\- n — h\ 

+ (1 - [cosh 20(q)][cosh 2e{q - k)]) ■ 

sr-^ nsMqa, q)] + nB[uJbiqo±, q - k)] ^ 
, , ^ - ^b{qo, q) - uJbiqo±, q-k) 
± [sinh 2^(q)][sinh 2^(q - k)] ■ 
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nsMgo, q)] + riBluJbiqoT, q - k)] \ 

„ -ujb{qo,q) -(^biqoT,q- k) j' 

(50) 

When computing goi above, we use the multiphcation table ±± = + and =i=± = — . The 
property Ub^qo, q) = UJb{qo, — q) displayed by the spectrum for Schwinger bosons ( l22ll 
was exploited to obtain the reduced expression (150|) . Expression (150|) can be used to 
show that {2n)~^ duj{S{-\-) ■ 'S"(+))|^;^ = at fc = 0. This is consistent with a spin- 
singlet hidden-order antiferromagnetic state inside the subspace where ^ • J2a ^z,i,a = 0: 
m2iJ2a '^hal'^) ~ ^- Expression fl50l) for the dynamical spin correlator is easily 
evaluated in the zero-temperature limit, where ideal BEC of the Schwinger bosons into 
the {qo,q) = (—,0) ground state occurs (126|) . At large sq, in particular, the third and 
last term of this expression is given by 

z{SH ■ S'{-))\ti = n-'A,{-, k){[ujb{-, k) - uj]-' + k) + CO]-'), 

for the case of hidden magnetic order, where the poles in frequency have spectral weight 

A^i-, k) = 7r(l - x)2(so/2)n2[tanh 2e{0)][Q^{k)/ui^\k)], 

with uil^\k) = [f^y(fe) — Q,'j_{k)]'/'^ . It reveals a Goldstone mode at = 0. Ideal 
BEC yields poles in frequency that exhibit a spin gap in the hidden channel (— ) at 
fc = from the first and second terms of expression fISUl) . on the other hand. The 
present mean-field theory for the hidden half-metal state therefore generally shows long- 
range antiferromagnetic order across the d-\- and d— orbitals, which is consistent with 
Goldstone's theorem. 

At large electron spin sq, the true self-consistent mean-field theory dynamically 
suppresses inter-orbital hopping, however. This leads to Schwinger bosons with 
degenerate spectra: uJb{qo,k) = uil^\k) for go = +,~- Expression (150|) in conjunction 
with ideal BEC then yields the result [T8] 

z{S{±) ■ S'{±))\i^ = 7r''A{±,k){[col'\k)-cor + U'\k)+u:]-'), (51) 

where the poles in frequency have spectral weight 

A{±, k) = 7r(l - xfsoh^[n^{k)/n±{k)]'/\ (52) 

Here, Q± = Q\\ ±Q±. The above dynamical spin correlator coincides with the transverse 
spin susceptibility, xj_(A;,a;), in the present zero-temperature limit by the fluctuation- 
dissipation theorem. Notice that the nature of the spectral weights fl52l) above implies 
the existence of hidden (— ) spin waves near zero 2D momentum with a spectrum ( H2i) . as 
well as true (-I-) spin waves near cSDW momenta with a spectrum (jH]) [18]. We conclude 
that the hidden half-metal state shows strict long-range antiferromagnetic order across 
the d-\- and d— orbitals of the iron atoms in the presence of inter-orbital hopping. 

The true spin excitations (+) predicted by the dynamical spin susceptibility flSTl) 
near the quantum-critical point are shown graphically by Fig. [3] for a set of parameters 
that is applicable to iron-pnictide high-temperature superconductors. Notice the spin 
gap (l45l) that exists at cSDW wavenumbers kcSDW = (7r/a)a; and {7i/a)y. It collapses 
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to zero at a critical Hund's Rule coupling — Jqc of moderate strength (H7|) . This strongly 
suggests that the hidden half metal state gives way to a cSDW metal phase [33] that 
shows strict long-range cSDW order at the QCP that separates the two phases. The 
critical cSDW metal has a low ordered moment because of proximity to the hidden half 
metal state. 



3.3. Fermi Surfaces 

We shall now obtain the electronic structure of the hidden half metal state by computing 
the one-electron propagator directly from the Schwinger-boson-slave-fermion mean- 
field theory for the two-orbital t-J model. It is given by the convolution of the 
propagator for Schwinger bosons with the propagator for slave fermions: iG{k,u) = 
Ea=d-,d+<^S° * Gi{^*|fc,^, where iG'^^l{ri,t]rj,t') = {fi^a{t)f]^p{t')). A standard 
summation of Matsubara frequencies yields the expression |18] 

^^nB[ub{q)] + npiSfiq - k)] 



fco=0,7r " ^ 



-UJb{q)+ef{q-k) 



+ 



r,nB[uJbiq)] + nplejiq - k)] 



+ [sinh e{q)Y ' ' . (53) 

^ '^'^ u + U!,iq)+ef{q-k) J 

Above, npis) = (exp[(;U — £)/kBT] + 1)^^ is the time-reversed Fermi-Dirac distribution. 
All of the Schwinger bosons condense into the lowest-energy state at 3 momentum 
(Q'o, q) = (tt, 0) as T — 7- by the principal meanfield equation (j26l) . This results in the 
following coherent contribution to the electronic spectral function at zero temperature 
and at large Sq: 

ImGcoh(fc, w) = soTi" ^ 5[uj + £f{ko,k)\. (54) 

fco=0,7r 

In the case of the unoptimized generic mean-field theory, it reveals inner and outer 
hole Fermi surfaces centered at zero 2D momentum that are depicted by Fig. O 
These collapse into doubly-degenerate circular hole pockets in the self-consistent mean- 
field theory at large sq, where inter-orbital hopping is dynamically suppressed. The 
contribution due to the Fermi-Dirac terms in the expression above represent incoherent 
excitations in the electronic structure. At energies u below the electronic Fermi level, 
the pole in the second term of expression fl53|) above represents the combination of a hole 
excitation, Sf > fi, with a spinwave, Ub > 0. This composite excitation therefore shows a 
gap ( 145|) at cSDW momenta, Acsdw- Notice that copies of the previous inner and outer 
Fermi surfaces now centered at cSDW wavenumbers {Ti/a)x and {TT/a)y exist at the 
QCP because of the collapse of the spin gap there: Acsdw 0. (See Fig. [3l) This is a 
spin-density wave nesting mechanism in reverse, where low-energy spinwaves centered 
at cSDW momenta produce nested Fermi surfaces! The incoherent contribution due to 
the first term in expression ( l53i) above vanishes in the zero-temperature limit at energies 
u below the electronic Fermi level. Last, the ratio of the incoherent spectral function 
integrated over momentum in the vicinity fc = or {-k /a)x{y) compared to the coherent 



Fermi- Surface Nesting in Reverse from Local-Moment Description 



17 



(A) IMAGINARY G (k.co) near QGP: ky = 



(B) IMAGINARY G (k,to) near QGP: ky = 
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Figure 4. Shown is tlie imaginary part of expression ()53|) with the parameter 
set that is listed in the previous caption to Fig. [31 but at Hund's Rule coupling 
—Jo ~ —Joc — O.OOIJ]^. The electron energy hu! is measured with respect to the top 
of the hole band, — e||(0). 



counterpart is Ylk^^^ ^i^cik, uj) / ^^lmGcoh{k,uj) = {-f/inso) ■ [{Au - A)Q\\/ (vo/a)^]. 
Here, Ao; measures how far in energy the hole hes below the Fermi level, A = in 
the region of 2D momentum centered at fc = 0, and A = Acsdw in the region of 2D 
momentum centered at kcSDW- It is assumed that Ao; > A. Also, 7 is the anisotropy 
parameter of the spinwave dispersion in question. (See fig. [3^.) 

We will now evaluate the former incoherent contribution to the spectral function in 
the large-So limit in the vicinity of the QCP, where A^sdw is small, at energies just below 
the Fermi level. As predicted by the mean- field theory, we shall assume that inter-orbital 
hopping is dynamically suppressed: ef{kQ, k) = e\\{k). The hole Fermi surfaces centered 
at zero 2D momentum (Fig. [2]) then reduce to two degenerate circular hole pockets at 
low doping. The previous long-wavelength approximations for the spinwave dispersion 
near zero 2D momentum fH2|) and near cSDW momenta fj44l) are then valid. Also valid 
is the long- wavelength approximation for the dispersion of the slave fermions (134 ah . with 
t\\ < 0. The imaginary part of the poles in the first (— ) and second (-I-) terms of Eq. 
( l53l) enforce energy conservation: —00 = ±ul^\q) + e\\{q — k). Taking the product of 
the two poles and enforcing energy conservation yields = [00 + e\\{q — fc)]^ — [a;^°^(q)]^. 
Substituting in the long-wavelength approximations for the spectra of the Schwinger 
bosons and of the slave fermions then yields a quartic equation for |q|. Its roots can be 
obtained analytically. We thus computed the incoherent (Fermi-Dirac) contribution to 
the electronic structure lmGinc{k,uj) in the thermodynamic limit from expression (l53ll 
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Figure 5. Shown is the imaginary part of expression (|53|) near the QCP at fixed 
energy just below the Fermi level, with the parameter set that is listed in the previous 
caption to Fig. H) 



by integrating the (5-function that enforces energy conservation over radial momentum 
\q\ analytically, and by performing the remaining angular integral numerically. 

Figure S] compiles the net electronic structure that is predicted by the present 
Schwinger-boson-slave-fermion mean-field theory. It is predominantly hole-type about 
zero 2D momentum because of the coherent contribution ( 15^ . which is independent of 
the Fermi level. Fig. also displays an electron-type "V" that separates purple from 
black, however. It meets the coherent hole band at the Fermi level, and its slope coincides 
with the longitudinal spin-wave velocity Vq. Figures andHt exhibit the energy gap 
below the Fermi level, A^sdw^ that exists at wavenumbers near {7T/a)x. A mix of 
electron and hole structure exists here at electron energies below the gap. The slope of 
the electron-type "V" that separates pink from purple in Fig. |3)d also coincides with 
the longitudinal spin- wave velocity vq, whereas the slope of the electron-type "V" that 
separates purple from black in Fig. coincides with the transverse spin-wave velocity 
vq/JcSdw- We conclude that the group velocity of the electron-type dispersion at 2D 
momentum near a point on one of the Fermi surfaces centered at a cSDW wavenumber 
(vr/a, 0) or (0,7r/a) is equal to the spin-wave velocity at the former momentum. Figure 
Hb also exhibits hole dispersion that follows the band u = —6\\[k — {7T/a)x], however. Its 
origin is the divergence of the coherence factor sinh^ 9{q) in Eq. fj53|) as the spin- wave 
frequency u^^\q) begins to vanish at g = {'K/a)x{y) near the QCP. Figure Ht shows 
that this effect is much less prominent along the transverse scan in 2D momentum. Last, 
Fig. [5] displays the same electronic structure as Fig. HI but at a constant energy that lies 



Fermi- Surface Nesting in Reverse from Local-Moment Description 



19 



halfway between the Fermi level and the minimum of the electron-type "V" dispersion 
at sSDW momenta. Panel (B) shows that traces of the hole-type dispersion centered at 
these momenta exist along the longitudinal {k^) axes. 

4. Exact Diagonalization and Conjecture 

The previous Schwinger-boson-slave-fermion meanfield theory for the the two-orbital t- 
J model ([1]) predicts that inter-orbital hopping is dynamically suppressed in the hidden 
half- metal state at large electron spin sq. Does this result persist for true electron spin, 
So = 1/2? Below, we will see that it does indeed persist when inter-orbital hopping 
is purely along next-nearest neighbors, while coherent inter-orbital hole propagation 
is possible in the opposite case where inter-orbital hopping is purely along nearest 
neighbors. 

We have obtained the exact low-energy spectrum of the two-orbital t-J model ([I]) 
for 4x4x2 local spin-1/2 moments plus one mobile hole by computer calculation. Each 
quantum state is specified by the combination of a spin background over all of the 4x4x2 
sites, confined to the subspace with total spin 5*2 = 0, and one hole location at down- 
spin sites in the spin background. The Heisenberg-exchange and Hund-exchange terms 
in the t-J model Hamiltonian ([1]) reduce to permutations of the spin backgrounds, and 
these are stored in memory. The matrix elements for correlated hopping terms in ([1]) are 
computed directly at each application of the Hamiltonian operator, on the other hand. 
Periodic boundary conditions are imposed, and translation and reflection symmetries 
on the 4x4x2 lattice are exploited in order to reduce the dimension of the Hilbert 
space. Global swap of the orbitals, P^j, is included in the list of symmetry operations 
at the extreme where inter-orbital next-nearest neighbor hopping is suppressed: = 0. 
For example, exploiting reflections about both principal axes of the 4x4 square lattice 
of iron atoms brings down the dimension of the Hilbert space to 75,624,211 states 
in such case at zero-momentum and even reflection parities. The low-energy states 
of the resulting block- diagonal Hamiltonian are obtained by employing the Lanczos 
technique |34j. We use the ARPACK subroutine library for this purpose [35]. in which 
case the Hamiltonian operation = Hlip) is accelerated by exploiting parallel threads 
with OpenMP directives. 

Figures El [7] and E\ display the evolution of the low-energy spectrum of one mobile 
hole in a 4 X 4 lattice of spin-1 iron atoms with the strength of the Hund's Rule coupling, 
—Jo- We chose the following set of Heisenberg exchange coupling constants and hopping 
matrix elements: j| = 0, > 0, j| = O.SJf*- = J^, tf = — 5Jj^, tj-{x) = —2J^, 
tiiv) = +2Ji", and t| = = t^. Recall that global swap of the orbitals, P^j, is then 
an exact symmetry. (See the end of section 2.) Also, the mean- field result ( H71) then 
predicts a quantum critical point at Jo = —0.8J^ in the thermodynamic limit {x = 0) 
that separates a hidden half-metal state at weak Hund's Rule coupling from a cSDW 
at strong Hund's Rule coupling. Figure [6] displays the exact spectrum in the absence of 
Hund's Rule. Red and blue points denote states that are respectively even and odd under 
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MOBILE HOLE IN HIDDEN FERROMAGNET WITH NO HUND COUPLING & tj^ = 
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Figure 6. Shown is tlie low-energy spectrum of tlie two-orbital t-J model, Eq. ([T]), over 
a 4 X 4 X 2 lattice with one hole, in the absence of Hund's Rule: Jq — 0. The remaining 
parameters are jj' 0, > 0, j| = 0.3 = J^-, if -?>Jt, (*) -2Ji^, 
tiiy) = +2J^, and 4 = = i^. A comparison with the hole spectrum, e/(fc) = 
and with the spin- wave spectrum, LOb{k) = ujl^\k), at large sq and a; = is also 
shown. Henceforth, some points on the spectrum are artificially moved slightly off 
their quantized values along the momentum axis for the sake of clarity. 



P^j. The spectrum is consistent with a hidden half metal. In particular, the lowest- 
energy states at fixed momentum carry spin-1/2 and they are nearly doubly degenerate. 
Furthermore, their dispersion agrees with that predicted by the previous mean-field 
theory (!22|) for the low-energy spin-wave excitations of the hidden half-metal state. 
Finally, the spin-1/2 states in Fig. [6] are well separated from the next excited states at 
fixed momentum, which carry spin-3/2. This observation is again consistent with well- 
defined spinwave excitations like those predicted by the dynamical spin susceptibihty 
( 15T|) of the Schwinger-boson-slave-fermion mean-field theory. 

Figure [7] shows how the spin-1/2 state at cSDW wavenumber (vr/a)^ with odd 
parity under P^j, in the absence of Hund's Rule, comes down in energy to become 
degenerate with the zero-momentum doubly-degenerate groundstate at a critical Hund's 
Rule coupling of Jqc = —1.73J^. We have measured the expectation value for local 
orbital swap at the hole iron site, Prfj(o). It has eigenvalues -|-1 and —1, which correspond 
to a hole with Sd^z and with 3dyz orbital character, respectively. Table [1] lists the 
corresponding expectation value over the groundstate at cSDW wavenumber {7i/a)x, 
(Prfj(o)) = —0.47, which means that the hole has 74% Sdy.^ orbital character. We 
note that turning off inter-orbital hopping entirely results in a somewhat higher critical 
Hund's Rule coupling of[l8] Jqc = —2.27J^. This dependence on inter-orbital hopping 
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Figure 7. The low-energy spectrum of the two-orbital t-J model with one hole roaming 
over a 4 X 4 X 2 lattice is displayed at the QCP that separates the cSDW from the 
hidden half-metal state: Jqc — —l-733J^. The remaining set of parameters are listed 
in the previous caption to the Fig. |6l 



should be compared to the mean- field result (H7|) . which predicts no dependence at all. 
The dispersion of the low-energy spin-1/2 states shown in Fig. [7]resembles the mean-field 
prediction for the dispersion of critical spin- wave excitations shown in Fig. [3l The former 
states are also well separated from the next excited state at fixed momentum. This again 
is consistent with the mean-field prediction ( 15T]) of well-defined spin-wave excitations at 
the QCP. Finally, the ordered moment for antiferromagnetic order across the d+ and 
d— orbitals that characterizes the hidden half metal and for the cSDW is defined by 
= [2/xb/(A/'fc-|)] ELo e*(^o"+'' '-')Si,„, with respective 3-momenta k = (vr, 0, 0) 
and k = (0, vr/a, 0). Table[T]lists the auto-correlation of each over the groundstate at zero 
2D momentum: {^i{k) ■ /c))o. They are given in units of the ordered moment of the 
ferromagnetic state, k = (0, 0, 0), over the 4x4x2 lattice: = (33/31)(2/is)^. Notice 
that ficSDW remains small at the QCP, at zero 2D momentum, which is consistent with 
neutron diffraction studies in iron-pnictide systems[7]. Notice also that fihFM remains 
sizeable at the quantum critical point, at zero 2D momentum, which indicates that 
hidden half metal character persists there. 

Last, Fig. [8] displays the exact low-energy spectrum of a true 4x4 lattice of spin- 
1 iron atoms with one mobile hole. In particular, Hund's Rule is enforced by setting 
Jo = — 23Jf'-. As before, red points and blue points are even and odd under global 
orbital swap, P^jj. We measured the expectation values for local orbital swap at the 
hole iron site, -Pdj(o), and we found that the hole in even and odd parity states has 
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Table 1. Listed are groundstate expectation values of physical observables for one 
hole hopping over a 4 x 4 x 2 lattice with the set of t-J model parameters used in Figs. 
[5][51 t2 = 0. Below, the QCP occurs at Jqc = — 1.733Jf'~, and the integer coordinates 
(ux, Uy) specify the momentum of the groundstate in units of 7r/2a. 



Observable Jq = @ (0, 0) QCP @ (0, 0) QCP @ (2, 0) Jo = -23 Jj^ @ (1, 0) 
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Figure 8. The low-energy spectrum of the two-orbital t-J model with one mobile 
hole hopping over a 4 x 4 x 2 lattice is displayed at strong Hund's Rule coupling: 
Jo — — 23Jf'". The remaining parameters are listed in the previous caption to Fig. [51 



orbital character that is respectively over 95% ?>dxz and ?>dyz. (See Table [TJ) We have 
hence replaced the former labels with the latter ones in the legend to Fig. El The 
groundstate notably has spin-1/2, it carries momentum {7r/2a)x{y), and it has orbital 
3dyz{3dxz) character. This is analogous to the groundstate momentum of {n/2a){x±y) 
that is predicted for one mobile hole in a 2D Neel state by Kane, Lee and Read [28]. 
Table [1] lists ordered moments computed in the groundstate at momentum {7T/2a)x{y). 
These moments combined with the low-energy spectrum indicate that the groundstate 
at thermodynamic hole densities x is a robust cSDW metal with Fermi surfaces that 
are centered at wavenumbers {'K/2a)x and {TT/2a)y. This state is therefore unable to 
account for any of the Fermi surfaces that are observed experimentally in iron-pnictide 
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Table 2. Listed are the characteristic properties of observable versus hidden quantum- 
critical spin-wave excitations in the two-orbital Heisenberg model that corresponds to 
the two-orbital t-J model ^ in the absence of mobile holes. (See ref. fTT|, Fig. 5.) 



spin wave 


order parameter 


2D momentum 


Hund's Rule? 


Pdd & 


observable 




(7r/a,0) & (0,7r/a) 


obeyed 


even 


hidden 


Si.d+ ~ Si,d- 


(0,0) 


violated 


odd 



FERMI SURFACES AT QCP & t/ = 




Figure 9. Conjectured nesting of Fermi surface pockets at QCP with nearest-neighbor 
inter-orbital hopping: and = 0. (Cf. figs. [2]and[7l) 



systems. 

The exact spectrum Fig. [7] at the QCP that separates a cSDW metal from a 
hidden half metal therefore indicates that coherent inter-orbital hole propagation exists 
near cSDW momenta in the case where inter-orbital hopping is purely across nearest 
neighbors. This contrasts with the prediction of dynamical suppression of inter-orbital 
hopping by the previous Schwinger-boson-slave-fermion mean- field theory at large sq- 
General agreement between the two calculations nevertheless exists. In particular, both 
the present exact results (Fig. [7j) and the previous Schwinger-boson-slave-fermion mean- 
field theory analysis find a QCP at moderate Hund's Rule coupling, —Jq, where one-hole 
groundstates at zero 2D momentum and at cSDW momenta become degenerate. The 
mean-field theory indicates that the cSDW momentum of the degenerate groundstate 
in Fig. [7] is carried entirely by a spinwave that softens to zero excitation energy at 
the QCP. (See Figs. El H] and El) Table El summarizes the nature of quantum-critical 
spinwave excitations for the corresponding Heisenberg model in the absence of mobile 
holesjH]. The key point to notice is that quantum-critical spinwaves at cSDW momenta 
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Figure 10. Shown is the low-energy spectrum of the two-orbital t-J model, Eq. ([T]), 
over a 4 X 4 X 2 lattice with one hole, in the absence of Hund's Rule: Jq ~ 0. The 



remaining parameters are J| 



0, > 0, j| 



0.3Jf 



Jn-L, t 



0, 



i!] = 0, t^ix + y) 



iJ^, and i; + y) — +iJ^- A comparison with the hole 



spectrum, Sf{k) = £||(fc), and with the spin-wave spectrum, ujb{k) 
So and a: = is also shown. 



,(0) 



(fc), at large 



are observable, and they hence have even parity under P^j, while those at zero 2D 
momentum are hidden, and they hence have odd parity under P^j. [See Eqs. fl^ and 
( 152|) . and Fig. |3l] We have confirmed this by exact diagonahzation of the corresponding 
Heisenberg model. The exact results displayed by Fig. [7] find one state at cSDW 
wavenumber [7r/a)x{y) that is degenerate with groundstate at zero momentum. It is 
odd (even) under P^j. Combining the meanfield theory picture of the QCP with the 
present exact results then leads to the following conjecture: the critical cSDW spinwave 
at wavenumber {ii/a)x relates the odd parity 3dyz portion of the hole pockets centered 
at zero 2D momentum (Fig. |2]) with its counterpart centered at the cSDW momentum, 
while the other critical cSDW spinwave at wavenumber {'ir/a)y relates the even parity 
3dxz portion of the hole pockets centered at zero 2D momentum with its counterpart 
centered at the other cSDW momentum. Further, we conjecture that the 3dxz orbital is 
related to the 3dyz orbital at the four hot spots in 2D momentum at which they coincide 
on the hole Fermi surface pockets by the hidden spin wave with odd parity under P^j at 
zero 2D momentum. Figure [9] summarizes the conjectured reverse nesting mechanism. 

Figures HD] and E] show the exact spectrum of the two-orbital t-J model (II]), 
respectively, in the absence of Hund's Rule coupling and at the QCP, but in the opposing 
case where inter-orbital hopping is purely across next-nearest neighbors. The model 
parameters, in particular, are j| = 0, Jf*- > 0, jj = O.SJf*- = J^, t\ = — SJf*-, t^ = 0, 
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Figure 11. The low-energy spectrum of the two-orbital t-J model with one hole 
roaming over a 4 x 4 x 2 lattice is displayed at the QCP that separates the cSDW 
from the hidden half-metal state: Joc = —2.26J^. The remaining set of parameters 
are listed in the previous caption to the Fig. 1101 



tl = 0, + y) = —iJi, and a; + y) = +iJi. Pink and light blue colors 

represent states that are respectively even and odd under the symmetry operation P^j 
that denotes swap of the d± orbitals after the gauge transformation e^*'^/^. The critical 
Hund's Rule coupling at the QCP is — Jqc = 2.26Jf'-, which is very close to its value 
in the absence of inter-orbital hopping[T8], — Jqc = 2.27 J^. Notice that both values 
are almost three times larger than the mean- field prediction at x = 0, — Jqc = O-SJ^. 
The quantum-critical spectrum shown by Fig. [11] is also very close to the corresponding 
one in the absence of inter-orbital hopping up to a rigid energy shift that is relatively 
small [18]. Further, the moments for hidden ferromagnetic order and for cSDW order at 
the QCP that are listed in Table [3] match those obtained previously in the absence of 
inter-orbital hoppingp^ to within 1%. We also computed the groundstate expectation 
values of modified orbital swap at the iron hole site, P^^{°), and these are listed in Table 
[3l A hole in a 3dx>z orbital has even parity (+1) under it, while a hole in a Sdy/z orbital 
has odd parity (—1) under it. Here, x' = (x + y)/2^^'^ and y' = (y — x)l2}l'^ are the 
2D coordinates along the next-nearest neighbor links. Notice that is generally 

small compared to unity, which means that the hole does not possess well-defined ?)dx'z 
or 'idyiz orbital character at the QCP. To conclude, good agreement exists between exact 
results for the spectrum of the hidden half metal in the presence of purely next-nearest 
neighbor inter-orbital hopping and dynamical suppression of the latter as predicted by 
the previous Schwinger-boson-slave-fermion mean- field theory at large sq. It suggests 
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Table 3. Listed arc groundstate expectation values of physical observables for one 
hole hopping over a 4 x 4 x 2 lattice with the set of t-J model parameters used in 
Figs. [TU]and[TTJ = 0. Below, the QCP occurs at Jqc = — 2.26Jf'", and the integer 
coordinates {nx,ny) specify the momentum of the groundstate in units of 7r/2a. 



Observable Jq = @ (0, 0) QCP @ (0, 0) QCP @ (2, 0) 
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Figure 12. Predicted nesting of Fermi surface pockets at QCP with next nearest- 
neighbor inter-orbital hopping: — Q and ^ 0. (Cf. figs. [^andlTT]) 

the reverse Fermi surface nesting mechanism that is outhned by Fig. [121 Last, it is 
useful to point out that related behavior is predicted theoretically for a mobile hole in 
the Neel state over the square lattice, where nearest-neighbor hopping (ti) is dynamically 
suppressed at low energy, leaving effective next nearest-neighbor hopping (t2) within the 
same antiferromagnetic sub lattice [3T| [36] [37]. 

5. Discussion and Conclusions 

Starting from a local-moment description of a cSDW over a square lattice of spin-1 iron 
atoms with mobile holes, we have succeeded in accounting for the nested Fermi surfaces 
pockets centered at zero 2D momentum and at cSDW momenta that are characteristic 



Fermi- Surface Nesting in Reverse from Local-Moment Description 



27 



of iron-pnictide high-temperature superconductors. In particular, zero-energy spin-wave 
excitations at cSDW momenta combine with hole Fermi surface pockets centered at 
zero 2D momentum to produce Fermi surface pockets centered at cSDW momenta. The 
former hole pockets exist because of proximity to a hidden half metal state with opposing 
polarized spin over 3d(^x+iy)z and 3d(^x-iy)z orbitals, respectively, which violates Hunds 
Rule. The isotropic {d-\-, d—) orbital basis that we choose notably maximizes the Hund's 
Rule coupling (see Appendix), and it leads to isotropic Heisenberg exchange coupling 
constants across neighboring spins on the square lattice of iron atoms. This orbital basis 
then very likely minimizes the net magnetic energy in the two-orbital t-J model ([1]) at 
fixed iron moment, (| Si^a]"^)^^"^ . Notice that the reverse Fermi-surface nesting effect 
uncovered here is a mirror image of weak-interaction descriptions for electronic structure 
in iron-pnictide high-temperature superconductors, where low-energy spin excitations 
centered at cSDW momenta are a result of Fermi surface nesting [2] |38j. Unlike current 
descriptions of cSDW's in iron-pnictide materials that start from the weak-interaction 
limit p] |38] . however, the present local-moment description predicts that hidden low- 
energy spin excitations across the d+ and orbitals exist near zero 2D momentum 
because of proximity to the hidden half metal state. 

It is useful to contrast our results with those of more ad hoc theoretical models 
that separate local moments from itinerant electrons [39] [IQ]. Although such models are 
capable of simultaneously accounting for the spin- wave spectra and for the Fermi surfaces 
seen in iron-pnictide high-temperature superconductors and their parent compounds, 
they clearly have less predictive power by virtue of the explicit separation between the 
two phenomena. Early models for iron-pnictide systems that simply add Heisenberg 
exchange interactions and Hund coupling to one-electron hopping Hamiltonians that 
already include nested Fermi surfaces, but that do not project out double occupancy 
at an iron site-orbital, also suffer from this drawback in our opinion|;4l]. In particular, 
such models essentially operate in the weak-interaction limit, but they fail to link Fermi 
surface nesting to low-energy spin excitations at cSDW momenta in iron-pnictide high- 
temperature superconductors. 

In summary, a mean-field theory analysis and an exact diagonalization study 
indicate that the two-orbital t-J model ([T]) for iron-pnictide high-temperature 
superconductors transits from a cSDW to a hidden half metal state with decreasing 
Hund's Rule coupling if off-diagonal magnetic frustration exists: e.g., J^^''^='= = Q, 
jd±,dT ^ g ^ showu by Fig. [H intra-orbital hopping stabilizes 

antiferromagnetic order across the d+ and d— orbitals, which characterizes the hidden 
half metal. This resembles the iron-pnictide phase diagram, where the cSDW is carrier 
poor while the superconductor is carrier rich. We therefore propose the hidden half 
metal state for the normal state of iron-pnictide high-temperature superconductors. This 
proposal is consistent with the low-energy spin-excitation spectrum and with the low- 
energy electronic structure shown by iron-pnictide high-temperature superconductors. 
Schwinger-boson-slave-fermion mean-field theory ( 15T|) predicts low-energy spinwaves 
that disperse anisotropically at cSDW momenta, (vr/a, 0) and (0, vr/a), near the quantum 
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critical point that separates the hidden half metal from the cSDW. The predicted 
dispersion of the spinwave spectrum, Fig. |3l notably shows a local maximum at the 
Neel wavenumber (vr/a, vr/a), which agrees with inelastic neutron scattering studies 
of the parent compound CaFe2As2 [12]. Further, the anisotropic dispersion that we 
predict at the QCP for low-energy spinwaves at cSDW momenta is consistent with 
recent observations of the same in iron-pnictide superconductors [Ti] [15] |16j |17j . 

Both the Schwinger-boson-slave-fermion mean-field theory formulation ( 153|) and 
exact diagonalization of the two-orbital t-J model ([T]) predict Fermi surface pockets 
around zero 2D momentum and around cSDW momenta at the quantum critical point 
if the intra-orbital hopping is hole type (see Figs. Hli, |5]and[7j). The Fermi surface 
pockets centered at (0, 0) are due to the coherent motion of holes in the hidden half 
metal state (see Fig. [2l). On the other hand, the low-energy electronic excitations about 
cSDW momenta are due to the combination of a hole excitation near the Fermi surfaces 
centered at zero 2D momentum with a low-energy spinwave that carries momentum near 
(vr/a, 0) or (0,7r/a). Figures Hb andHb show that the electronic structure is a mixture 
of electron and hole types near those points in 2D momentum space. The dispersion 
velocity of the electron-type feature is notably given by the spin-wave velocity there! 
The above predictions are not inconsistent with the electronic structure shown by iron- 
pnictide high-temperature superconductors (cf. refs. [5] and |6|). The general features 
of the electronic structure that we predict near the quantum-critical point that separates 
the cSDW from the hidden half metal with the bare minimum of two 3d orbitals are 
also similar to what is predicted by band-structure calculations that include all five 3d 
orbitals. In particular, the checkerboard structure of the inequivalent iron atoms in 
iron-pnictide materials results in zone-folded Fermi surface pockets at (7r/a,7r/a) that 
have only weak spectral weight [38] . 
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Here we calculate the Hund's Rule coupling in the 2D subspace spanned by the 3dxz and 
Sdyz orbitals in iron-pnictide materials. Consider the most general pair of basis states 
for such orbitals: 

<p{r) = i?3,2(r)[(cosao)e-^'^«r2,+i(^^) + (sin ao)e+^'^°r2,-i(^^)], (0.1a) 
^(r) = R,,2{r)[-{smao)e-"^'Y2,+,{n) + (cosao)e+^'^«F2,-i(^^)]. (O.U) 

Notice that ao = or 7r/2 corresponds to the 3d(^x+iy)z/'id(^x-iy)z orbital basis, whereas 
Qo = 7r/4 and 0o = corresponds to the '2>dxzl^dyz orbital basis. The exchange Coulomb 
integral is related to the Hund's Rule exchange coupling constant Jo by 



dVi / d^r2<p*{ri)i){ri)- 
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in general. In the present case, this yields 
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where VL' is the solid angle rotated by 00 about the z axis. The integrals over solid 
angles and ^2 can be performed in the standard way|42] by use of the mathematical 
identity 
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Performing the remaining radial integrals after substitution of the hydrogenic radial 
wave function 
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are the strengths of the integrals due to the / = 4 and / = 2 channels, respectively. Note 
that the binomial-like series that appear above result from the difference 

(5-/)! -2-(^-')4'r'^(l/2), 



X 



where Sn (x) is the m*'^ derivative of the finite geometric series sum Sn{x) = ^^=o 
The Hund's Rule coupling is then 

1 

— Jn = (cos"^ an + sin'' an). (0.9) 

30.3082 2ao^ u u; k j 

Notice that cos^ ao + sin'^ ao = 1 — | sin^ 2ao, which reaches its maximum value of unity 
at ao = or 7r/2 and its minimum value of 1/2 at oq = it /A. We conclude that the 
Hund's Rule couphng is largest in the isotropic 3d(^x+iy)z/'id(x-iy)z orbital basis, while it 
is smallest in the anisotropic ^dxz/^dyz orbital basis from Chemistry. 

It is useful to compare the maximum Hund's Rule coupling in the isotropic 
X — iy)z orbital basis with the on-site Coulomb repulsion in that case: 

f/o = j dh^j ^^242(^)1^2.1(1^1)! ^^ '1^^^ 42(^2)1^2,1(^^2)1^. (0.10) 

Addition of angular momentum yields the identity 

\Y2,im^ = ^PoicosO) + -^P2(cos^^) - ^P^{cose) (0.11) 
in zoTT in 

in terms of Legendre polynomials. Substituting it above, with P;(cos6') = [An/ {21 + 

1)] ^^^^,0(^)5 along with the mathematical identity (10. 3p . yields the result 
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for the Coulomb integral, where 



12' /25 2^ 2^ 2^ 2 1 \ 
= i± ^ _ 5^ + 10^ - 10- + 5- - - (0.13) 
" 2i2\^ 7 8 9 10 11 127 ^ ' 

is the strength of the integral in the / = channel. The ratio of the on-site Coulomb 
repulsion fl0.12p to the Hund's Rule coupling (10. 6p is then 

- — ^ = 10.6743. (0.14) 
-'0 

Study of the exchange integral (10. 2p yields that this ratio coincides with the ratio 
jf"'"''^'*'(drct)/jf'''''^~(drct) in the regime ao ^ a. 
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